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Abstract 

For a reliable prediction of the NMSSM Higgs boson signatures at present and future high-energy 
colliders and a proper distinction of the NMSSM and MSSM Higgs sector the precise knowledge 
of the Higgs boson masses including higher-order corrections is indispensable. In this paper, the 
one-loop corrections to the neutral NMSSM Higgs boson masses and mixings are calculated in 
three different renormalisation schemes. In addition to the DR renormalisation scheme, existing 
in the literature, two other schemes are adopted. Furthermore, the dependence on the value of 
the top quark mass is investigated. The resulting Higgs mass corrections have been compared 
and the residual theory error due to missing higher-order corrections can be estimated to be of 
the order of 10%. 
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1 Introduction 



Supersymmetric theories [Tl[2] provide a natural solution to the hierarchy problem [2l[3]. The latter 
is related to the fact that the Standard Model (SM) Higgs boson mass receives radiative corrections 
which are quadratic in the cut-off scale, usually chosen to be the GUT scale Mqut = 10^^ GeV. 
To keep the Higgs mass of the order of the electroweak (EW) scale an extreme fine-tuning of 
the model parameters is necessary. Supersymmetry (SUSY) introduces a new symmetry between 
fermionic and bosonic degrees of freedom, leading to extra contributions to the radiative corrections 
of the Higgs boson mass due to the new SUSY particles which cancel the dangerous quadratic 
divergence of their corresponding SM counterpart. In this way the Higgs mass is naturally kept at 
a phenomenologically valid level even in the presence of high mass scales. 

As the superpotential must be analytic in the chiral superfields, two complex Higgs doublets 
have to be introduced, Hu to provide masses to the up-type fermions and Hd to ensure non-zero 
down-type fermion masses. In this way the theory is also kept anomaly-free. In the superpotential 
of the Minimal Supersymmetric Extension of the Standard Model (MSSM) [1] the two Higgs doublet 
fields mix through the term fiHueH^ which involves the higgsino mass parameter fi. Phenomenology 
requires fi to be of the order of the electroweak scale [5]. On the other hand, /j, is present in the 
superpotential before the symmetry breaking and hence not a priori of the order of the EW scale. 
In the Next-to-Minimal Supersymmetric Extension of the Standard Model (NMSSM) [BHS] this fi 
problem is solved by dynamically relating the value of /i to the EW scale. The parameter arises 
here as the vacuum expectation value of the neutral component of an additional Higgs field S which 
is a singlet field with respect to the SM gauge groups. It couples to the MSSM Higgs fields via 
the interaction term \S{Hu€Hd). With the scalar field S acquiring a non-zero vacuum expectation 
value (VEV) Vs/V2 an effective term is generated, naturally expected to be of the order of the 
EW scale, /i = \{S) = Xvg/V^- Furthermore, new contributions to the quartic coupling increase 
the tree-level mass value of the lightest Higgs boson. In the MSSM the tree-level mass of the 
lightest Higgs boson is predicted to be below the mass of the Z boson. Large radiative corrections 
involving top and stop loops are necessary to lift the Higgs mass value beyond the lower bound 
from direct searches at LEP [10]. The additional NMSSM contributions lift the mass of the SM-like 
Higgs boson more easily beyond the LEP bound. 

The extension of the Higgs sector by two more degrees of freedom through the introduction 
of an additional singlet superfield field ^ leads to a total of 7 Higgs bosons after electroweak 
symmetry breaking, three neutral CP-even, two neutral CP-odd and two charged Higgs bosons. 
The fermionic component of S mixes with the MSSM higgsinos and neutral gauginos to yield five 
neutralinos. Among the considerable phenomenological modifications in the Higgs and neutralino 
sector compared to the MSSM are possible new Higgs-to-Higgs decays, as e.g. the decay of a SM- 
like scalar Higgs boson into a pair of lighter pseudoscalar Higgs states, so that the present Tevatron 
and LHC search studies for supersymmetric Higgs particles have to be revisited [11]. Such a Higgs 
boson could have escaped the LEP bounds [12] . 

From the above discussions it is clear that the precise knowledge of the Higgs boson masses is 
indispensable to distinguish between MSSM and NMSSM predictions, to properly define scenarios 

^We denote superfields by a hat over the field. Fields without hat are the corresponding component fields. 
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with new Higgs-to-Higgs decays within the NMSSM and to correctly interpret the experimental 
data. In the MSSM, where radiative corrections are crucial to accommodate the light Higgs boson 
mass with the LEP limits, the Higgs boson masses have been calculated with impressive accuracy. 
After the calculation of the dominant one-loop contributions due to top/stop loops [13] the full 
one-loop corrections have been provided in |14j . The leading logarithmic two-loop effects obtained 
through renormalisation group equations [15] have been completed by the genuine 0{atas) [T6H20] . 

[MllinillT], 0{ahas) [22] and 0{atai, + a^) [23] two-loop corrections in the limit of zero 
external momentum. The electroweak two-loop effects including part of the external momentum 
dependence have been calculated [24] and the leading three-loop contributions have been evaluated 
in [25]. Also in the case of a CP-violating MSSM a great effort has been undertaken to calculate 
the higher-order corrections. After first investigations [26], they have been evaluated at one-loop 
order in the effective potential approach [27] and with the renormalisation group improved effective 
potential method through next-to- leading order [28]. In the Feynman diagrammatic approach 
the one- loop leading mf corrections have been provided in [2^, a full one- loop calculation in 
j30] . and more recently the leading two-loop contributions of 0{atas) have been evaluated in 
|31j . The corrections have been implemented in the public computer code FeynHiggs [18 1 130 1 
I32j . CPsuperH [33], another public code, is based on the renormalisation group improved effective 
potential approach [28]. The MSSM Higgs mass spectrum for real parameters can also be obtained 
from the spectrum calculators |34fl36j . 

In the NMSSM, however, the higher-order calculations to the Higgs boson masses have not yet 
reached the same level of accuracy as in the MSSM. The leading one-loop contributions due to 
top/stop and bottom/sbottom loops have been calculated in the effective potential approach [37], 
the one-loop contributions due to chargino, neutralino and scalar loops have been evaluated in 
leading logarithmic order in Ref . [38] . These corrections and the leading logarithmic two- loop terms 
of 0{atas) and 0(af), taken over from the MSSM results, have been implemented in the public 
computer code NMHDECAY [39]. Furthermore, the full one-loop contributions have been computed 
in the DR renormalisation scheme and the 0{atas + OibOs) corrections have been provided in the 
approximation of zero external momentum [30]. They have been implemented in NMHDECAY as well 
as in the spectrum calculator SPheno [36]. 

In this paper we complement the effort to reach a higher level of accuracy in the computation of 
the NMSSM Higgs boson masses by providing the full one-loop corrections in a mixed DR— on-shell 
and in an on-shell (OS) renormalisation scheme. By comparison with the results of Ref. [40] in 
the DR scheme, the dependence on the renormalisation schemes can be studied. In this way an 
estimate of the theoretical error due to missing higher-order corrections can be derived. Having 
available the corrections in the mixed and in the OS renormalisation scheme in addition enables a 
comparison of the NMSSM results in the MSSM limit with the corresponding MSSM corrections 
given in the OS scheme. 

The paper is organised as follows. In section [5] we present the details of our calculation starting 
by introducing the NMSSM Higgs sector and setting up our notation in 12.11 After presenting the 
chargino and neutralino sector in 12.21 we list the parameters, which we employ, in section 12.31 The 
renormalisation is explained in detail in 12. 4[ The explicit computation of the one-loop corrected 
Higgs boson masses and mixing matrix elements is described in 12. 5[ Section [3] finally is devoted to 
the numerical analysis. Our results are summarised in[31 
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2 Calculation 



2.1 The NMSSM Higgs Boson Sector 

The Higgs mass matrix is obtained from the NMSSM Higgs potential, which is derived from the 
NMSSM superpotential Wnmssm, the corresponding soft SUSY breaking terms and the D term 
contributions. The NMSSM superpotential for the Higgs superfields Hu,Hd,S in our conventions 
reads 

Wnmssm = Wmssm - eabXSH^Ht + ^^tS^ , (1) 

where a,b = 1,2 are the SU{2)i fundamental representation indices and eab denotes the totally 
antisymmetric tensor with ei2 = e^^ = 1. The Higgs superfield which couples to up- type (down- 
type) fermion superfields is given by -ffu(d)- The parameters A and k are dimensionless and, working 
in the CP-invariant NMSSM, are chosen to be real. The first term in Eq. ([T|) is the MSSM 
superpotential, which in terms of the quark and lepton superfields Q, U^, D'^, L, read^ 

Wmssm = eabiVeHSt'E' + VdHSQ^D' - VuKQ^'U'] , (2) 

where we have omitted colour and generation indices. Note that the MSSM ^ term in NMSSM 
constructions is commonly assumed to be zero as well as terms linear and quadratic in S. The soft 
SUSY breaking terms in the NMSSM in terms of the component fields Hu,Hii, S are given by 

Csoft = Csoft,MSSM - + {eabXAxSH^Ht - ^kA^S^ + h.c.) , (3) 

with the MSSM soft SUSY breaking Lagrangian 

^soft,MSSM = -m\^H\Hd-m\^HlHu-mQQ'^Q -m\VL-mljU^^^ 

- mle*ReR - {e,b[yeAeHSL% + yaAaH2Q''d*R - y„A,i?^Q^n^] + h.c.) 

- ^{MiBB + M2WiWi + M3GG + h.c.) , (4) 

where in the first two lines tilde denotes the scalar component of the corresponding quark and 
lepton superfield, and Q = {uL,dL)'^, L = {i'L,eL)'^- Note that e.g. u*j^ is the scalar component 
of U^. In the last line the soft SUSY breaking gaugino mass terms for the gaugino fields B, Wi 
{i = 1,2,3) and G are given. All soft SUSY breaking trilinear couplings Ak {k = X,K,d,u,e) and 
gaugino mass parameters Mj (i = 1,2,3) are assumed to be real, and squark and slepton mixing 
between the generations is neglected. Furthermore, possible soft SUSY breaking terms linear and 
quadratic in the singlet field S are set to zero in accordance with the majority of phenomenological 
NMSSM constructionfil. 

The neutral components of the Higgs fields can be parametrised in terms of CP-even and CP-odd 
fiuctuations around their vacuum expectation values, 

= -^{vd + hd + ittd) , Hi = ^{vu + hu + iau) , S = -^{vs + hs + ius) , (5) 



^The superscript c denotes charge conjugation. 

^In the MSSM there is an additional soft SUSY breaking term m^^jj^eabH^H^ with m^^H^^ usually expressed 
in terms of the soft SUSY breaking parameter B and the higgsino parameter /i as rnjj^fj^ = _B/j. This term is 
not explicitly added to the NMSSM soft SUSY breaking Lagrangian but is generated dynamically via the term 
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with the VEVs Vd,Vu,Vs chosen to be real and positive. The parameter /i then arises dynamically 
from the singlet field expanded about its VEV, cf. Eq. ([T]), 



V2 



(6) 



The minimisation conditions of the tree-level scalar potential can be applied to replace the soft 
SUSY breaking Higgs mass parameters ^^/f^ , "i^^ , m| by combinations of A, k, the electroweak 
gauge couplings g and g' , the VEVs and the trilinear couplings Ax, A^, so that the tree- level mass 
matrix M| of the neutral CP-even Higgs bosons obtained from the second derivative of the Higgs 
potential with respect to the fields in the vacuum, in the basis = {hd, h^, hg)'^ can be cast into 
the form 



Mi 



( 



(A^ - g^)vuVd 



}?VdVs - (Rx + RVs)Vu 

where we have defined 



V 



(A^ - f)VuVd - {Rx + RVs/2)Vs X^VdVs - {Rx + RVs)Vu\ 

X^VuVs - {Rx + RVs)Vd 

2t^vl + 



X^VuVs - {Rx + RVs)vd 



(7) 



+ RkVs j 
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l2\ 



R = \k . 



R> 



Ri^ — --j=kA^ 



(8) 



Note, that the MSSM limit can be recovered by A,k — )• (with the ratio k/X kept constant for a 
smooth approach) and keeping the parameter /i = Xvg/V^ as well as the parameters A\ and 
fixed. In this limit we hence have — )• 00. 



The CP-even mass eigenstates Hi 



1,2,3) are obtained by an orthogonal transformation 



7^ (the summation over paired indices is implicit). 

Hi = 'R'ijhj , .7 = d,u,s . 



(9) 



They are ordered by ascending mass with < Mh^ < Mh^- The CP-odd fields (0^,0^,0^) can 
be rotated first to separate a massless Goldstone boson G 





( 






^ ad 


as = 








I] 


au 


G I 








\ as 



(10) 



with the abbreviations = cos a; and Sx = s'mx, in the following also = tanx. Starting from 
the tree-level CP-odd mass matrix squared in the basis {ad,au,as)^ and using Eq. (jlOp results in 
the mass matrix squared Mp in the basis = {a, as, G)"^, which reads 



Mi 



{2Rx + Rvs)vs^ 



20 



{Rx - RVs)vCAl3 



{Rx - Rvs)vcA(5 {2R + Rx/vs, ^ 
\{Rx + Rvs/2)vs-^f {Rx - Rvs)vsAfs 



(Rx + Rvs/2)vs^f^ 
SRi^Vs {Rx - RVs)vSAl3 

{2Rx+RvsW^ I 



(11) 
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with A/S = P — /Sn and v'^ = v^^ + v^. The angle /3„ coincides at tree- level with the angle (3 defined 
through the ratio of the VEVs Vu,Vd, tan 13 = Vu/vd, Hence A/? = at tree-level leading to a 
massless Goldstone boson, which decouples as it should. The first entry {Mp)ii at tree-level, 

iMj,)n = + , (12) 

S2fS 

becomes the mass of the pseudoscalar Higgs boson in the MSSM hmifl Applying an orthogonal 
rotation TZ^ to = {a,as,G)'^ , the CP-odd mass eigenstates Ai = Ai,A2,G [i = 1,2,3) are 
obtained, 

Ai = nf^hf , (13) 

where at tree- level TZ^^ = 1 and TZ^^ = TZf-^ = for i 7^ 3. The pseudoscalar masses are ordered by 
ascending mass, Ma^ < Ma2- The CP-even and CP-odd Higgs mass values squared are given by 
the eigenvalues of the respective mass matrices Mj and Mp. Analytic expressions would be rather 
complicated if not expanded in special parameter regions as e.g. in Ref. [JT], where a comprehensive 
investigation of the NMSSM Higgs boson sector has been performed. In our analysis the mass 
eigenvalues are derived numerically. The charged Higgs boson mass on the other hand takes a 
simple form after the massless charged Goldstone boson has been separated by an orthogonal 
rotation with a mixing angle /3c. At Born level we have f3c = (3 and get 

Ml^ = M^, + V2A ^{A, + K^) - ^ = + (M|)n - ^ , (14) 

■52/3 V2 ^ 2 

with Mw being the W boson mass and where we have applied the definition given in Eq. (jl2p . 

2.2 The Neutralino and Chargino Sector 

The chargino sector remains unchanged with respect to the MSSM. For completeness and to set up 
our notation we briefly repeat the chargino system. The chargino mass matrix depends on the wino 
mass parameter M2, the effective Higgsino parameter /i = \vs/^/2 and the ratio of the vacuum 
expectation values tan /?. In the interaction eigenbasis it is given by 



Mc = ( r""' ^^^^^ ) . (15) 
It is diagonalised by two real matrices \J and V ^ 

WMcV^ U = 0. and V=[ ifdetM^O 

I (J3O+ if det Mc < ' ^ ' 

with the Pauli matrix CJ3 to render the chargino masses positive. The rotation matrices 0± are 
given in terms of the mixing angles 

, 2V2Mw{M2Cp + piSp) 2^Mw{M2S^ + iJiCp) 

tan2C7„ = T — r> — — , tan2y_i_ = ^ — ^ — — , 17 



A\ takes the role of the soft SUSY-breaking parameter B in the MSSM. 



6 



and the two chargino masses read 

^\t2 = I + + 2M^ T mi - i?f + AM^{M'^c% + M| + + 2M2/xs2/3)]5 } . (18) 



m 

xr 



In the neutrahno sector, the mixing of the fermionic component of the singlet superfields S 
with the neutral gauginos B,Wz and higgsinos H^,H^ yields in the Weyl spinor basis ip^ = 
{B, W3, H^, H^, S)'^ the 5x5 neutralino mass matrix 



Mn 



f Ml -cpswMz spswMz \ 

M2 cpcwMz -spcwMz 

-cpswMz cpcwMz -\vs/V2 -\vu/\/2 

spswMz -spcwMz -Xvs/V2 -Xvd/V^ 

\ -Xvu/V2 -Xvd/V2 V2kvs J 



(19) 



We have introduced the short-hand notation sw = sin^vK, cw = cos^vK for the Weinberg angle 
9w, and Mz denotes the Z boson mass. Diagonalisation with a unitary matrix M yields the five 
neutralino mass eigenstates Xi {i = •■•i^), 

X°=AA,,V'°, j = l,...,5. (20) 

They are ordered by ascending mass, m^o < ... < m^o. With our assumption of CP invariance and 
allowing for negative neutralino mass eigenvalues, the mass matrix N is real. 

2.3 Parameter Basis 

The NMSSM Higgs potential depends on 12 independent parameters in the CP-conserving case. 
They are given by the soft SUSY breaking mass parameters itl^h^ , 1 ) the gauge couplings 
g, g' , the vacuum expectation values Vu, Vd, Vg, the dimensionless parameters A, k and the soft SUSY 
breaking couplings AxjA^. For the physical interpretation it is convenient to replace some of these 
parameters: The minimisation of the Higgs potential V requires the terms linear in the Higgs fields 
to vanish in the vacuum. Hence for the scalar fields, 

where the brackets denote the vacuum. The corresponding coefficients, which are called tadpoles, 
therefore have to be zero. At lowest order the tadpole conditions for the CP-even fields reaqj 

\dh^/ ^ th^ = Vd-rriH^ - RxVuVs + Vd{vd - v J + ^Mvu + Vs) - ^VuV^ = (22) 

= thu=Vum%^-RxVdVs + ^-^^Vu{vl-vl) + YVu{vl + Vs)-^Vdv'i = 0{23) 

/ dv\ 



''As we work in the real NMSSM the derivatives of the Higgs potential with respect to the CP-odd fields are zero 
and no additional conditions have to be required. 
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Equations (f22l [23l can be exploited to replace it^'h^,^'h^,^'s by the tadpole parameters th^, t^^ 
and ths- The parameters g,g',Vu,Vd are replaced by the electric charge e, the gauge boson masses 
M\Y, Mz and the ratio of the vacuum expectation values tan/3 = tjj = Vu/vd through the relations 

eMz , eMz 

9 = I 9 = -jTj— (25) 

Jul - 



l + tl eMz "^""V l + tl eMz- 



Finally, Ax is replaced by the charged Higgs boson mass Mjj± , 



with 



2Mw L Ml 



w 



Ml ' 



(28) 



where we have kept the dependence on the tadpole parameters and the mixing angles f3c = Pn = (^B- 
At tree-level they coincide with (3. It is the angle (3 which will be needed for the renormalisation: 
According to our renormalisation schemes, presented in section [27^ the renormalisation procedure 
can be performed before the transformation into mass eigenstates. This means that the mixing 
matrices, also those separating the neutral and charged Goldstone bosons from the Higgs bosons 
and hence the angles /3„ and /3c, which appear in these matrices, do not receive counterterms. 
Therefore, special care has to be taken to determine the elements of the mass matrices without 
inserting the tree- level relation Pb = P and to apply the renormalisation procedure uniquely for 
tan /3 given by the ratio of Vu and Vd- 

To summarise, we work with the following parameter set 

thu, th„ th^, e, M^, M|, tan/3, M^±, A, k, v^, . (29) 



2.4 Renormalisation Schemes 

For the determination of the loop-corrected Higgs boson masses the Higgs self-energies have to be 
calculated. They develop ultraviolet (UV) divergences. Evaluating the self-energies in D = 4 — 2e 
dimensions, the divergences can be parametrised by 1/e leading to poles in D = 4 dimensions. 
To get a finite result the parameters entering the loop calculation have to be renormalised by 
the introduction of appropriate counterterms absorbing the UV divergences. Note, that at higher 
orders also the terms linear in the Higgs fields get loop contributions. Therefore, also the tadpole 
parameters th. (i = u, d, s) have to be renormalised, in order to fulfill the tadpole conditions 
Eq. (I2ip . Before we describe the determination of the one- loop corrected Higgs boson masses in 
section 12. 5| the renormalisation schemes which have been adopted, shall be presented in detail in 
the following. 

We start with a renormalisation scheme which is a mixture between on-shell and DR renormal- 
isation conditions. In order to make contact and to compare to earlier results presented in Ref. |4U| 
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our result will also be converted to a pure DR scheme. To get an estimate of the uncertainties 
due to missing higher-order corrections we study the influence of the renormalisation scheme on 
the one-loop Higgs mass corrections. We therefore also compare to a pure OS scheme. In all three 
renormalisation schemes, mixed, DR and OS, we start out from the parameter set Eq. (j29p and 
replace the parameters by the renormalised ones and their corresponding counterterms, 



Ml, 



M| + 5Ml 
M^ + 5M^ 
Ml, + 6Ml, 
{l + SZe)e 
X + S\ 

Vs + 5Vs 



thd tha + Stha 

tan /3 —7- tan (3 + 5 tan (3 

n n + 5n 

— >■ Aj^ + 5Afi . 



(30) 



For the field renormalisation, the Higgs boson doublet and singlet fields are replaced by the renor- 
malised ones and a corresponding single field renormalisation constant for each doublet and the 
singlet, respectively, 



(31) 



5 



ZsS 



1 + -5Zs ] S . 



Applying this renormalisation procedure, the renormalised self-energies in the basis {/if 



1,2,3} = {hd,hu,hs} for the CP-even Higgs bosons and in the basis {hf 
for the CP-odd ones can be derived as {X = S, P) 

1.1 7-7 /'-II It-' 



1,2,3} = {o,o„G} 



! 3 

1 



[5Z,.,^x{Ml)kj + {Ml),k6Z^xj^x] - 5{Ml)ij with i,j = 1,2,3 , (32) 



where 



and 



5Z, 



hfhf 



6Z, 



hlhl 



SZh^, 6Z^ 



hlhl 



SZs, dZ, 



for i / j. 



(33) 



^^hPhP 



sl6ZHa + c15Zh^ , 



^^hPhP 
^^hPhP 



spcp{6ZHa 



5Z, 



hPhP 



SZs , 



5Z, 



hPhP 



c}6Zh, + s}6Zh^ (34) 
(35) 
(36) 



Here, M| and MJ, are the tree-level CP-even and CP-odd Higgs boson mass matrices squared, 
respectively. For the derivation of the counterterm matrices (5(M|) and 5{Mp), the CP-even and 
CP-odd Higgs boson mass matrices squared are expressed in terms of the parameter set Eq. ([29]) 
including also their dependence on the tadpole parameters as well as on the mixing angle Pb- As 
the expressions are quite lengthy they are given in Appendix A. Then, the parameters entering 
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these matrices are replaced according to Eqs. (130]) and an expansion about the counterterms is 
performed. The part of and MJ,, respectively, which is linear in the counterterms, corresponds 
to 5(M|) and 5{Ml). 

The renormalised self energies in Eq. ()32p are related to the ones in the basis of the mass 
eigenstates through, cf. Eqs. (j9ll3p . 

W.(fc') = 7^ffc7^f^^:,|,s(fc2) (37) 

= 7^ffc7^j;%^p(A;2) with i,i,A;,/ = 1,2,3. (38) 
The field renormalisation constants 6Zhj^,SZh^,SZs are obtained from 

szh^h, = \nl^?6ZH, + \nU''5ZH^ + , i = i, 2, 3 , (39) 

with 



2\ div 



i = 1,2,3, (40) 

fc2 = (M(J'))2 



where (M^J',^)^ denotes the corresponding tree-level mass squared. The field renormalisation con- 
stants are defined in all three renormalisation schemes via DR conditions. This is indicated 
by the superscript 'div' and means that in the field renormalisation only the divergent part 
A = 2/(4 — D) — 7s + ln(47r) is kept with being the Euler constant. Solving Eq. (|39|) for 
5Zh^ , 5Zh^ , 5Zs results in 

5Zh, = [(I^f3l'l^f2l' - \nU^\T^U^)5ZH,H, + (|7^f2p|7^f3p - \nU^\nU^)5ZH,H, 

+ (IT^fsl'l^fsl' - \nf^\''\nl^\'')5ZH,H,]/Rf (41) 
5Zh^ = [(|7efi|2|7ef3|2 - |7^f3|2|7^f,|2)^z^^^^ + {\nU''\nl,\'' - \nf,\''\nU'')5ZH,H, 

+ (l^fll'l^fsP - \Rf^?\nl^?)5ZH,H,]/R'r (42) 

5Zs = [(|7^f2p|7^fl|2 - |7^fl|2|7^f2|2)5z^,,^,, + {\'Rf^?\Rl^? - \nU^\nl^\^)5ZH,H, 

+ (l^fsl'l^fiP - \Rf^\''\nl^\'')5ZH,H,]/Rf (43) 

where 

jyS IX**^ I^IT?*^ I^IT?*^ 1^ _1_ IX**^ I^IX**^ I^IX**^ 1^ _1_ IX**^ IRIX'S' l^lx**^ 1^ 

rt^ — I '"'-111 l''t'22l I '■'-33 1 "T I '■''23 1 I '■'-32 1 "r|'"'-12l I ''^21 1 I '''-33 1 

- I^f2pl^f3l'l^fll' - I^f3l'l^fll'l^f2l' + l^f3l'l^f2l'l^flP • (44) 



Contrary to the field renormalisation constants the renormalisation conditions for the remaining 
parameters are different in the three chosen renormalisation schemes, as will be described in the 
following. 

Mixed renormalisation scheme 

In the mixed renormalisation scheme we divide the parameters into parameters defined through 
on-shell condition J§ and into parameters defined via DR conditions: 

Mz, Mw,MH±,th^,th^,th,,e, tan /3,X,Vs,K,A^ . (45) 

' ' ' V ' 

on-shell scheme DR scheme 

^In slight abuse of the language we also call the renormalisation conditions for the tadpole parameters on-shell. 
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In the following, the various counterterms shall be specified in more detail. 

(ijii) Gauge boson masses 

The gauge boson masses are defined through on-shell conditions, 

ReSi^(M|) = 0, Ketl^wiM^) = , (46) 

where T denotes the transverse part of the respective self-energy. For the mass counterterms this 
yields 

<5M| = ReS^^(M|), SM^, = ReY.^w{M^) . (47) 

Note, that the NMSSM gauge boson self-energies differ from the MSSM case due to the introduction 
of an additional superfield S. 

(in) Mass of the charged Higgs boson 

The mass of the charged Higgs boson is determined through the on-shell condition, 

ReEH±/fT(M|±) = 0, (48) 

resulting in the corresponding counterterm 

(5M^± = ReS H± (M2 ± ) . (49) 

(iv-vi) Tadpole parameters 

The tadpole coefficients are required to vanish also at one-loop order, yielding 

4? - ^^h, = ' i = d,u,s , (50) 

where t^/^^ stands for the contributions coming from the corresponding genuine Higgs boson tad- 
pole graphs. As the tadpole graphs are calculated in the mass eigenstate basis, they have to be 
transformed to the interaction basis. Applying Eq. ([9]) we have 

6th,=n^it'^H] > i = d,u,s, j = 1,2,3. (51) 



(vii) Electric charge 

The electric charge is defined to be the full electron-positron photon coupling for on-shell ex- 
ternal particles in the Thomson limit, so that all corrections to this vertex vanish on-shell and for 
zero momentum transfer. The counterterm for the electric charge is then given in terms of the 
transverse part of the photon-photon and photon-Z self-energies |43|^ . 



6Z, 



1 dJ^'^^jk') 



+ ^^^. (52) 



^Note that the sign of the second term in Eq. (|52|) differs from the one in [43] due to our conventions in the 
Feynman rules. 
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(via) tan /3 

For the renormalisation of tan /? we adopt the DR scheme. Applying Eq. (j3ip and 

Vi ^ Vi + 5vi i = u,d (53) 

results in 

tan/3 



S tan /9 = tan f5 











Vd ) _ 


div 



(54) 



div 



where we have used in the last step dvu/vu\div = Svd/vd\div |44] . The field renormalisation constants 
SZh^ and 6Zh^ are given in Eqs. ([ill Ii2|) . 

(ix) Coupling A 

In the mixed renormalisation scheme A is defined as a DR parameter. The counterterm is 
determined via the renormalised self-energy S^p^^p, see Eq. (j32]) . using that 



div 



5(M 



p)ii 



'hPhP 



(M 



div 



(55) 



As (5(Mp) 11 contains the counterterm Eq. ([55]) can be solved for 6\ resulting in 



6X 



4AM2,>^2^ 



div 



(56) 



The self-energy Tij^pj^p is obtained from the self-energies in the mass eigenstate basis S^.^. (i, j 

11 '''3 

1, 2, 3) through 



(57) 



(x) Singlet vacuum expectation value Vg 

The vacuum expectation value Vg of the singlet field is renormalised in the DR scheme. The 
counterterm is derived by exploiting the chargino sector. In fact, the lower right entry of the 
chargino mass matrix Mq in the interaction basis reads, cf. Eq. (jl5p . 



(MC)22 = ^ 



(58) 



Applying Eq. ()30p . expanding around the counterterms and extracting the terms linear in the 
counterterms yields 



— S{Mc)22-Vs^ 



(59) 



div 



'Compare with Eq. ([12)1 . 
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with (5 A given by Eq. ()56p . The counterterm d{Mc)22 is obtained from the renormahsed chargino 
self-energies in the following way. Defining the general structure of a fermionic self-energ}l§ as, 

where Vifi = (iT 75)/2 are the left- and right-handed projectors, we use the condition 



div 







(61) 



Note that is a 2 x 2 matrix with the entries given by the renormahsed self-energies of 

the charginos in the mass eigenbasis. The renormahsed self-energies in terms of the unrenor- 
malised ones, the mass counterterms and field renormalisation constants are given in Appendix B, 
Eqs. (|122p - (|125p . The structure of the condition Eq. (|6ip has been chosen such that the (divergent) 
contributions of the chargino field renormalisation constants drop out. Replacing the renormahsed 
self-energy by the relations (jl22p - (jl25p leads to the counterterm S{Mc)22, 



6{Mc)22 = \ 



(Mc)22[V^^Sj±(fc2)y + C/^Sj±(A;2)C/*]22 + [^7^Sji(A:2)y + ytsJ|(A;2)C/*]22]^,^(62) 

Note that the k? dependence in Eq. (I62p drops out as only the divergent part is taken. 
(xi) Coupling k 

The counterterm for the DR renormahsed parameter k is derived from the neutralino sector 
in an analogous procedure as the chargino sector was exploited to determine 5vs- The lower right 
entry {Mn)55 of the neutralino mass matrix Eq. (fT9]) reads 



(M;v)55 = \/2, 



leading to the counterterm 6k, 



6k 



-^^6{Mn)55 - K— 

V2Vs Vs 



(63) 



(64) 



The determination of 6vs has been described in the previous paragraph. For the determination of 
(5(MAr) 55 we use the condition 



(M^)55[AA^(S^o(fc2) + S«o(fc2))AA]55 + [M^{t^J{k'') + S«f (A;2))AA]55 



div 



0, 



(65) 



where the fermionic self-energy structure Eq. (I60p has been applied for the decomposition of the 
renormahsed neutralino 5x5 self-energy matrix S^o. Once again the condition has been chosen 
such that the divergent parts of the field renormalisation constants cancel in Eq. (j65p . Rewriting 
the equation in terms of the unrenormalised self-energies, c/. Eqs. (jl26p - (jl29p . yields 



<5(M^)55 = ^ ^ 



(M;v)55[AA^(S^o(A;') + S^o(fc2))AA]55 + [AA^(S^o^(A:=^) + S^^(A:=^))AA]55 ,. , (66) 



which is inserted in Eq. 



to determine 6k. 



div 



The decomposition can be applied both for unrenormahsed self-energies E and renormahsed self-energies E. 
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(xii) Trilinear coupling 

The trilinear coupling A,^ is also defined as a DR parameter. For the derivation of the coun- 
terterm we use that 



'hPhP \ iMp)22 



. (67) 



div 

Equation (j67|) depends on SA^ via the mass matrix squared counterterm 6{Mp)22- Solving for SA^ 
we have 



SA^ 



^ [W((M|)22) -5/] -M^ + ^-?]].. ■ (68) 



3kVs ' "2 "2 ^'^V "J '^^ K V. 

The counterterm 5f is derived from 



div 



J -z;?n2 + ^h^^p^PB\ + ^2„2^2 — y^H± - MwC^p\ 



+ ^MMkfll [2XM^^s'^s,p + ?,^e\l] , (69) 



with A/3 = 13 — Pb- This is done by replacing the parameters th^, th^, th^, e, Mz, M\y, M^±, tan/3, 
A, K, and Vg according to Eq. (j30p . performing an expansion about the counterterms, extracting 
the part linear in the counterterms and finally applying the tree-level relations for the tadpole 
parameters t^^ = t^^ = t^^ = and for the mixing angle, Pb = P- 

The self-energy T,f^Pf^p in terms of the corresponding self-energies in the mass eigenbasis is given 

by 

^h^hp = nf^^A^A.n^ , i, j = 1, 2, 3 . (70) 

Alternatively, we could have derived the counterterms 5vs and dn from the Higgs sector instead 
of resorting to the chargino and neutralino sector. We have explicitly verified that this leads to 
the same results for the one-loop corrected Higgs boson masses. Our choice of renormalisation 
allows a non-trivial cross-check of the renormalisation procedure. Moreover, it paves the way for 
an extension of the one- loop corrections to the Higgs decays into charginos and neutralinos. 

On-shell renormalisation scheme 

In this renormalisation scheme we keep the conditions (i)-(viii) (tan /3 is still renormalised in 
the DR scheme) but the parameters Vs,X,k and Af^ are determined via on-shell renormalisation 
conditions. This is done by applying on-shell renormalisation conditions on the two CP-odd Higgs 
boson mass eigenstates ^1,^21 on the mass eigenstates of the two charginos and on the mass 
eigenstates of the two lightest neutralinos. As OS conditions are imposed on the mass eigenstates 
not only single elements of the counterterm mass matrices occur in the conditions, in contrast 
to the DR conditions exploited within the mixed scheme. Thus, the equations get more involved 
and include also the parameters Mi,M2. This is why six renormalisation conditions are required, 
although in the end the counterterms for Mi, M2 are not needed. 

Requiring on-shell masses for the CP-odd Higgs bosons as well as for the two charginos and for 
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the lightest and next-to-hghtest neutrahnos leads to the following relations, 



(7^^ 5{Ml) (7^^)T 

V* 5{Mc) W 
V* 5{Mc) 
[M 5{Mn) 
[m 5{Mn) 



11 



11 



ReS^,^,((il/4°J)2), 



22 

(m^) + Sj± (m^)) + S^i (m^) + Sj| (m^) 



1, 



22 2 

_ 1 

11 ~ 2 

_ 1 

22 ~ 2 



Re 
Re 
Re 



. X2 ^ X22 X2' X22^ X2" X22^ X2' X22^ X2 ' 



\ All Xl All Al ' All Al All Al 



m 



X2^ X22 X2'^ Xoo^ Xo^/ Xoo^ Xo^ Xo, ^ Xo-^ 



X22 ' X2 



X22' X2' 



X22' xv 



(71) 

where (Mj^*^ )2 are the tree-level masses squared of the pseudoscalar mass eigenstates A\^Ai. The 
counterterm matrix 5(Mp) is derived as described above. The entries 6{Mc)ij = 1,2) of the 
counterterm matrix 6{Mc) for the charginos read 



<5Mc,i 
6Mc,2 

6Mc2, 
6MC22 



6M2 
6M\ 



8/3 



w 



+ (5tan/3 \/2c|Mvk 
(5tan/3 V2(?pMwsp 



V2 V2 



(72) 
(73) 

(74) 
(75) 



And finally, the entries 5{Mi\f)ij = 6{MN)ji = 1...5) of the counterterm matrix 6{Mn) for the 
neutrahnos can be cast into the form 



6Mn^^ =5Mi 



2Mzsw 



+ (5tan/3 



C[sMzSwS2l5 



5MjVi4 = (<5M| - 6M^) 
6MN22 =5M2 



+ otanp ■ 



SMn2, 

6Mn2, 



=6M, 



2 C/3 



w 



2Mw 



2Mzsw 
(5tan/3 c^Mvi^s^ 



cw 



6Mi 



2 «/3 



W 



2Mw 
V2 \/2 



(5tan/3 c^Mi^ 



(76) 
(77) 

(78) 

(79) 
(80) 

(81) 
(82) 



^/2eM]vsw 



(5tan/? 



V2XciMwsw 



-+ 



+ {X5Ze - 5X) 



V2MwswSi3 



5M, 



N45 



V2eMwsw 



V2e 



+ {X5Ze - 5X) 



\/2cj3Mwsw 



(83) 



(84) 
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JMat^s =V26kvs + V2k6vs (85) 
6Mn,2 =SMn,.^ = 6Mn2s = SMns-^ = SMn^^ = . (86) 

With these ingredients the exphcit expressions for the equations ()7ip can be derived. In our 
calculation, the system of equations (j7ip is solved numerically, keeping, however, the dependence 
on the divergent part A explicitly. 

PR renormalisation scheme 

The DR renormalisation scheme differs in the conditions (i)-(vii) from the mixed one, which 
means in the conditions for Mz , Mw , Mff± , th^ , i/i^ ,thg,e. Note that a change of the renormalisation 
condition for Mjj± can be interpreted as a change of the condition for Ax, cf. Eq. ([H]). For these 
parameters instead of OS conditions DR renormalisation conditions are adopted now. The other 
renormalisation conditions do not change. 



2.5 Loop Corrected Higgs Boson Masses and Mixing Matrix Elements 

The one-loop corrected scalar Higgs boson masses squared are extracted numerically as the zeroes 
of the determinant of the two-point vertex functions , 



±H,HAk') ^H,H,{k') - (m(?])2 + ±H,H.Ak') J 

In the same way the pseudoscalar masses squared are obtained from , 



(87) 



k^-{Mflf + tA,AAk^) tA,A,{k^) 

tA.AAk^) k^-{Mf>f + tA,AAk^) 



The superscript (0) denotes the tree-level values of the masses squared. It should be noted that in 
Eq. ()88p the mixing with the Goldstone bosons is not taken into account. We have checked explicitly 
that the numerical effect is negligible. The unrenormalised self-energy and tadpole contributions 
that occur implicitly in Eqs. (f87|) and ([88]) are evaluated at one-loop order. They contain fermion, 
Goldstone and Higgs boson, gauge boson and ghost loops as well as loops from the corresponding 
superpartners i.e. sfermions, charginos and neutralinos. 

The mass eigenvalues are obtained iteratively. In order to obtain the lightest scalar Higgs boson 
mass e.g., in the first iteration the external momentum squared k"^ in the renormalised self-energies 
^HiHj is set equal to the lightest scalar tree-level mass squared. Then, the mass matrix part of T"^, 
meaning {iT^ + /c^l), is diagonalised and the resulting mass eigenvalues squared are used in the 
next iteration where k'^ is set equal to the lightest of the obtained mass eigenvalues. Once again the 
mass eigenvalues are obtained. The procedure is repeated until the deviation between the lightest 
eigenvalue and the one of the previous iteration is less than 10~^. The other Higgs mass eigenvalues 
are derived accordingly 

Due to the radiative corrections, not only the masses of the particles receive contributions but 
also the fields are affected. To take these effects into account, new matrices, TZ^'^\ TZ^'^\ are 
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introduced which transform the fields hu,hd,hs and a, as into the corresponding one-foop mass 
eigenstates, respectively. These matrices are no physical observables and beyond lowest order 
they depend on the external momentum in the self-energies. For the derivation of the radiatively 
corrected matrices, 1Z^'^\ we follow the procedure applied in Ref. [30]. It ensures the correct 

on-shell properties for the external particle in processes with external on-shell Higgs bosons at 
higher orders and thus accounts also for the mixing between the Higgs bosons. This leads to finite 
wave function correction factors. In the scalar case e.g. we have to apply the additional factor Zj to 
the tree-level matrix TZ^ , which rotates the interaction states {hd,hu,hs)^ to the mass eigenstates 
{Hi,H2,H3f, to get the one-loop matrix elements, 

7^,^'^' = (Zf),,■7^f; , i,j = HuH2,H3, l = hd,hu,hs. (89) 

The correction factor is given by 

(Zf),, = , (90) 

with 

Zi = . (91) 

l + (ReSf)'(MHj2 

The prime denotes the derivative with respect to k'^, and M^, is the one- loop corrected Higgs boson 
mass squared. The effective self-energy T,^^ appears in the diagonal Higgs boson propagators 

It is given by (no summation over 

The off-diagonal Higgs boson propagator Aij [i ^ j, no summation over reads 

"p5 "p5 "p5 ■pS' 

A--(p) = iUi ZLJ! (OA) 

; psrs rs j_r)f^SrSrS _rSr2S _rs ■pS'i _rSrS2 ' ^ ' 

^ ii^ jj^ II ^ '^^ ij^ jr li ii^ jl jj li U ij 

The argument k"^ in F^- has been dropped for better readability. For Zfj we have in terms of the 
propagators (again no summation over the indices) 



Aij{k 



2\ 



Au{k^) 

£.,(M^J [mi - (Mg))^ + tujMf,^)) - MMh)MM%:) 

- [mI^ - {mP^? + %AMl)) [mi - {mPj^ + tu{Ml)) 
Zi = \. (96) 
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p 11 

In case of the pseudoscalar 2x2 mixing matrix TZ-j^' the effective self-energy reduces to 




(97) 



and the off-diagonal propagator reads 



a jj ij 



(98) 



The thus derived mixing matrix elements include the full momentum dependence and imaginary 
parts of the Higgs boson self-energies. Alternatively we could have set k'^ = which corresponds 
to the result in the effective potential approximation and yields a unitary mixing matrix. For our 
parameter sets used in the numerical analysis we found that the differences in the two approaches 
are negligible. Furthermore, the imaginary parts of the mixing matrix elements are small compared 
to the real parts. 

3 Numerical Analysis 

The calculation of the Higgs and gauge boson self-energies, of the tadpoles and the counterterms 
has been performed numerically in two different calculations. In the first calculation all necessary 
Feynman rules have been derived from the NMSSM Lagrangian and implemented in a FeynArts 
model file [55]. In the second calculation the Feynman rules have been obtained with the Mathe- 
matica package SARAH [36]. The Feynman rules in the two approaches have been cross-checked 
against each other and also against the rules given in Ref. [8]. Subsequently, in both calculations 
FormCalc |47j was used to evaluate the self-energy and tadpole diagrams in the 't Hooft-Feynman 
gauge, in which the Goldstone bosons and the ghost fields have the same masses as the correspond- 
ing gauge bosons. The divergent integrals are regularised applying the constrained differential 
renormalisation scheme |48j which has been shown to be equivalent [l9] to the SUSY conserving 
dimensional reduction scheme jSOj . The numerical computation of the integrals has been performed 
with LoopTools [JT]- Two Mathematica programs have been written to evaluate the counterterms, 
diagonalise numerically the one-loop corrected Higgs boson mass matrices and extract the mass 
eigenvalues. 

For our numerical analysis we follow the SUSY Les Houches Accord (SLHA) [51] and use as input 
values the Fermi constant Gp = 1.16637 • 10~^ GeV~^ and the Z boson mass Mz = 91.187 GeV. 
For the electroweak coupling we set a = 1/137. From these input values we derive the parameters 
of our input set defined in Eq. (j29|) . The top quark pole mass is given by Mf = 173.3 GeV. As we 
cross-check our results against the ones of Ref. [40] which uses the running DR quark masses, we 
need to calculate these as well. In order to obtain the DR top quark mass we convert Mf at the 
scale Qt = Mt in the corresponding running mass. The SM renormalisation group equations are 
then used to evolve the top mass up to a common scale Q chosen to be of the order of the SUSY 
breaking scale, where the gluino corrections are added. We denote the running top mass by mt in 
the following. The same procedure is applied to the bottom mass starting from the SLHA input 
value mf,(mfe)^^ set equal to 4.19 GeV. The masses of the light quarks are chosen as m„ = 2.5 MeV, 
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rric = 1.27 GeV, = 4.95 MeV, = 101 MeV [52]. The r mass has been set to nir = 1.777 GeV. 

In the fohowing we wih discuss the results for various scenarios which exemphfy different effects 
of the higher-order corrections. For our scenarios, we took care not to violate unitarity bounds 
by choosing X, k such that VA^ + < 0.7. Furthermore, Vg has been chosen to be of the order 
of the vacuum expectation value v. As the one-loop corrections to the pseudoscalar masses are 
small, we mostly show plots for the scalar masses and comment briefly on the pseudoscalar masses. 
Furthermore, the corrections for the heaviest scalar Higgs boson are not shown, as in all cases they 
are negligible. 

Variation of A We first verify that we reproduce the results of Ref. [10] by adopting a DR renor- 
malisation scheme and choosing the same parameter set 

K = A/5 , tan/3 = 2 , = 500 GeV , = - W GeV , = 250 GeV (99) 

with A being a free parameter and Vg given by Vg = ypiyujX. For the squarks and sleptons a common 
soft SUSY breaking mass Ms = 300 GeV has been adopted, and the remaining soft SUSY breaking 
parameters have been chosen as 

At = At, = Ar = -1.5Ms , Ml = Ms/3, Ms = 2/3M5 , M3 = 2Ms . (100) 

The renormalisation scale has been set equal to Qq = Ms and for the top and bottom quark mass the 
running DR mass at the scale Qq has been used. Furthermore, the light quark masses have been set 
to zero as in Ref. [IQ]. We find agreement for the one-loop corrected Higgs boson masses. Starting 
from this scenario in the DR scheme, in order to investigate the effect of different renormalisation 
schemes, the DR input values have been converted to the mixed renormalisation scheme as defined 
above, as well as to a pure on-shell scheme, and the Higgs boson masses and the matrix TZ^ are 
evaluated accordingly. For the definitions of the schemes, see section 12.41 In Fig. [T] we show the 
matrix elements squared {TZf^)^ and (7^23)^, respectively, of the mixing matrix 7^^. These matrix 




Figure 1: The matrix element squared {TZfs)^ (left), (T^fa)^ (right) for the two lightest CP-even Higgs 
bosons as function of A at tree-level (yellow/full) and at one-loop level adopting a DR (blue/dotted), a 
mixed (red/dashed) and an on-shell (green/small dotted) renormalisation scheme. 
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Figure 2: The mass Mhi of the hghtest (left) and of the next-to-hghtest (right) CP-even Higgs boson 
as function of A at tree- level (yellow/full) and at one- loop level adopting a DR (blue/dotted), a mixed 
(red/dashed) and an on-shell (green/small dotted) renormalisation scheme. 

elements are a measure of the strength of the singlet component of the two lightest Higgs bosons. 
They are shown for the three different renormalisation schemes compared to the tree-level result 
as a function of A. In order to match the result of Ref. [30] they have been evaluated at vanishing 
external momentum k"^ = 0. At tree-level, for small values of A, the lightest CP-even Higgs boson 
is dominantly MSSM-like and the next-to-lightest is dominantly singlet-like. With increasing A the 
mixing increases and at tree-level there is a cross-over at A ~ 0.94: the next-to-lightest Higgs boson 
is now more MSSM-like than the lightest Higgs boson. The higher-order corrections change the 
amount of the singlet component and for this parameter set there is no cross-over below A = 1. 
Furthermore, we see that for A > 0.2 the amount of the singlet component is hardly affected by the 
renormalisation scheme. The curves for the three renormalisation schemes lie on top of each other 
for A > 0.65. For smaller values of A, however, after having crossed (T^fs)^ = and {TZ^s)"^ — 1' 
respectively, (TZf^)^, ('T^fs)^ ™ OS-scheme start to differ largely from the corresponding values 
in the two other schemes. The reason is that the finite parts of the counterterms involved in the 
conversion from the DR parameters to on-shell parameters contain a division by A, so that these 
finite counterterm contributions blow up in the limit A — t- 0. The Hi mass squared even turns 
negative for A < 0.1. The matrix elements and masses in the OS scheme are therefore not plotted 
any more for A values below this value. 

Figure [2] shows the one-loop corrected mass Mhi of the lightest (left) and of the next- 
to- lightest (right) CP-even Higgs boson as a function of A in the three different renormalisation 
schemes compared to the tree- level result. The tree- level masses increase with rising A due to the 
NMSSM contribution ~ A^ sin^ 2/3 from the Higgs quartic coupling. As Hi is dominantly MSSM- 
like its one-loop corrections are much more important than for the singlet dominated H2. For the 
latter they are negligible at small A and more important for large A where H2 is more MSSM-like. 
As can be inferred from the figures the deviations between the different schemes is negligible. The 
curves for all three schemes lie on top of each other. Apart from small values of A where the one- 
loop corrected Hi and H2 mass in the on-shell scheme start to differ from the DR and the mixed 
scheme. 
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Figure 3: The matrix element squared (T^-fg)^ (left), (T^-fg)^ (right) for the two lightest CP-even Higgs bosons 
as function of A at tree-level (yellow/full), at one- loop level with the top quark pole mass (blue/dotted) and 
with the running DR top quark mass (red/dashed). 



The mass corrections and hence also the strengths of the singlet and MSSM component, re- 
spectively, strongly depend on the value of the top quark mass, which reflects the fact that the 
main part of the higher-order corrections stems from the top sector. Figure [3] shows the mixing 
matrix elements squared (7^f3)^, ("T^fs)^ as functions of A at one-loop level, calculated in the DR 
scheme, with the top quark mass taken as the running DR mass mt = 150.6 GeV in one case and 
as the pole mass Mt = 173.3 GeV in the other case. For comparison the tree-level values are shown 
as well. Figure m displays the corresponding Higgs masses. In contrast to the case where the top 
quark mass has been set to the DR value, using the top pole mass leads to a one-loop corrected 
lightest Higgs boson which is singlet-like and to a next-to-lightest Higgs boson being MSSM-like. 




Figure 4: The mass Mh^ of the lightest (left) and of the next-to-lightest (right) CP-even Higgs boson 
as function of A at tree-level (yellow/full), at one- loop level with the top quark pole mass (blue/dotted) and 
with the running DR top quark mass (red/dashed). 
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Figure 5: The matrix element squared {TZf^)^ (left), (T^-fa)^ (right) for the two lightest CP-even Higgs bosons 
as function of A at tree-level (yellow/full), at one-loop level at the renormalisation scale Qo = 150 GeV 
(blue/dotted), 300 GcV (red/dashed) and 600 GcV (green/smaU dotted). 

Furthermore, the Higgs mass corrections are more important for a higher top quark mass value. 
Defining the relative correction AMh /Mh as 

AMh \Mh-mP\ 

Mh mP ' ^ ^ 

where Mu {M^) denotes the 1-loop corrected (tree-level) Higgs boson mass, the relative correction 
for Hi amounts to maximally 55% for the running mass and 79% for the pole mass. The large 
corrections also explain the change of the one- loop corrected Hi from a MSSM-like to singlet-like 
Higgs boson. In fact, due to the large corrections the one- loop corrected Hi mass gets shifted above 
the one-loop corrected mass of H2. Due to our convention to label by ascending indices the Higgs 
bosons with increasing mass, the Hi one-loop corrected mass is assigned to H2 and vice versa, so 
that Hi and H2 interchange their roles and Hi becomes singlet-like at 1-loop whereas H2 becomes 
MSSM-like. 

In order to get an estimate of the missing higher-order corrections we investigate the influence 
of the renormalisation scale Qq . The results are shown (in the DR scheme) in Fig. [5] for the mixing 
matrix elements squared. They are plotted as a function of A at tree-level and at one-loop level for 
three different values Qo = 150, 300 and 600 GeV. Note that the scale Qq also changes the value of 
the running b and t quark masses. The corresponding plots for the masses of the two lightest Higgs 
bosons are depicted in Fig. [6l Whereas the change of the renormalisation scale alters the mixing 
matrix elements considerably, the effect on the Higgs boson masses is less pronounced. With rising 
Qq the running top quark mass decreases, leading to smaller one-loop Higgs boson mass corrections. 
The relative correction for the lightest Higgs boson mass is 63% for Qq = 150 GeV and 50% for 
Qq = 600 GeV. The relative correction of the next-to-lightest Higgs boson mass is changed less, 
with AMH2/MH2 = 15% and 8% for Qq = 150 and 600 GeV. This is due to the higher tree-level 
mass value. The residual theoretical uncertainties due to missing higher-order corrections can thus 
be estimated to 0(10%). 
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Figure 6: The mass Mhi of the hghtest (left) and of the next-to-hghtest (right) CP-even Higgs boson 
as function of A at tree- level (yellow/full), at one- loop level at the renormalisation scale Qo = 150 GeV 
(blue/dotted), 300 GeV (red/dashed) and 600 GeV (green/small dotted). 



The amount of singlet component of the Higgs bosons strongly affects their couplings to fermions 
and gauge bosons and hence their phenomenology. In particular, with the small Hi and H2 mass 
values between ~ 90 and ~ 180 GeV, the question arises if the scenario has been excluded by LEP, 
Tevatron or LHC. We have explicitly verified that there are still regions in A which have not been 
excluded. The allowed and excluded regions are shown in Fig. [71 The exclusion limit for A < 0.6 
where Mhi ^ 113 GeV is due to the LEP exclusioE0of in the channel e^e~ — t- ZH — > Zbb |53j . 
The Tevatron |54j and the present LHC results |55p56] do not constrain the scenario. Note, that our 

'^^The exclusion regions always apply to the one-loop corrected Higgs boson masses. 
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Figure 7: The Higgs boson masses at tree-level (dashed) and one- loop (full) for Hi (blue/dark grey) and H2 
(red/light grey) in the DR scheme with the exclusion limits set by LEP. 
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A« [GeV] A« [GeV] 

Figure 8: Left: The matrix element squared {TZfr^)^ [i ~ 1,2) as function of at tree-level (dashed) and 
one-loop (full) for Hi (blue/dark grey) and H2 (red/light grey). Right: The corresponding tree-level and 
one-loop corrected masses. 

limits represent a rough estimate and cannot replace a sophisticated study of exclusion limits set 
by a combination of the experimental results. We have cross-checked though our exclusion limits 
against those obtained with HiggsBounds [57j and have found agreement. 

Concerning the pseudoscalar Higgs bosons, the lighter state Ai is singlet-like, whereas the 
heavier one is MSSM-like, and the corrections are small. They are negative and in this scenario 
the relative corrections are below 2%. The masses increase with rising A. The one-loop corrected 
mass of the lighter pseudoscalar Higgs boson ranges between 40 and 140 GeV, and M^^ takes 
values of 586 to 598 GeV. 

Variation of In the following the variation of the soft SUSY breaking coupling parameter A^ 
is investigated. We have chosen a common soft SUSY breaking squark mass rriQ = 1.1 TeV and 
Ms = 600 GeV, large enough to fulfill the present exclusion limits on the squark masses of the first 
two generations and the gluino mass set by the LHC experiments [58j. The full parameter set is 
given by 

A = 0.6 , K = A/3 , tan /3 = 2 , = 500 GeV , ^ = 275 GeV 
At = Ab = Ar = -I.5M5 , Ml = Ms/3 , M2 = 2/3M5. , M3 = 2Ms . (102) 

The top and bottom quark masses have been chosen to be the pole masses, with mP°'° = 4.88 GeV. 
The renormalisation scheme is the mixed scheme as defined in section 12.41 The mixing matrix 
elements squared quantifying the amount of the singlet component of Hi and H2 are presented in 
Fig. [S] (left). They have been obtained by the procedure described in section [231 We have explicitly 
verified that the difference to the values obtained by setting fc^ = is negligible. The tree-level 
and one- loop corrected Hi and H2 masses are shown in Fig. [8] (right). As can be inferred from 
the figure, for large negative values of A,^ the lightest Higgs boson is mostly singlet-like and the 
heavier one MSSM-like. With increasing Af^ the mixing increases developing a rapid cross-over at 
A,^ w —294 GeV. The one- loop corrections shift the cross-over to a larger value Af^ w —240 GeV. 
This behaviour is also reflected in the Higgs boson masses. Below (above) the cross-over the 
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lightest (next-to-lightest) Higgs boson is mostly singlet-like and exhibits a strong dependence on 
A,^, whereas in the MSSM-like case the Higgs bosons hardly depend on Af^. This behaviour results 
from the fact that the soft SUSY breaking term which contains is a cubic coupling in the singlet 
field S. A variation of can hence only be communicated through a singlet contribution of the 
Higgs boson fields. Since the mixing is small, only fields which are mainly singlet-like are affected. 
Despite the singlet character of Hi for small A,^ the one-loop corrections can be important and 
more than triple the mass. This is mainly due to the small tree-level value. In case of the next-to- 
lightest Higgs boson H2 the one-loop corrections in the singlet-like case are less pronounced and of 
0(10 GeV). The corrections to the MSSM-like Higgs bosons increase the masses by about 37 GeV. 
No Higgs boson mass values Mh^ 2 have been excluded so far by the experiments. With a slightly 
higher value of A, however, the MSSM-like Higgs boson mass gets shifted above 141 GeV, so that 
the scenario would be excluded by the new LHC exclusion limits |56j . 

For the pseudoscalar masses the relative corrections are very small, at the 6% level for the 
singlet-like Ai and below 0.2% for the heavier A2. The one-loop masses decrease with decreasing 
absolute value of A^ and are 110 GeV < M^, < 345 GeV and 639 GeV < Ma^ < 640 GeV for 
GeV > > -400 GeV. 

Variation of Mjj^ The variation of Ax corresponds to a variation of the charged Higgs boson 
mass, cf. Eq. (jl4p . The other parameters are kept fixed and chosen as 

A =0.65, K = A/3, tan/3 = 2, A^ = -10 GeV , /i = 225 GeV 
At = Ab = Ar = -1.5Ms , Ml = Ms/3, M2 = 2/3Ms , M3 = 2Ms . (103) 

The common soft SUSY breaking squark mass is once again taken to be mo = 1.1 TeV and 
Ms = 600 GeV. The parameter A\ is varied between 250 and 600 GeV, which corresponds to 
Mjj± ~ 420. ..610 GeV. Outside this parameter range the Hi mass squared becomes neg ativeEl. 
The results for the singlet components and Higgs boson masses at tree-level and one-loop are 
shown in Fig. M The lightest Higgs boson Hi is dominantly MSSM-like with a 100% MSSM 
component at Mjj± ~ 520 GeV. Here H2 is maximally singlet-like, though not completely. The 
heaviest scalar Higgs boson takes over the remaining singlet component (not shown here), so 
that 3(^^13)'^ = 1 as demanded by unitarity of the mixing matri^lj • The one- loop corrections 

increase the mixing between Hi and H2 away from the maximum and minimum singlet values. They 
lead to a more pronounced maximum and minimum in the singlet component and slightly shift their 
positions to smaller H^ masses, Mff± ~ 518 GeV. The Hi and H2 masses are maximal and minimal, 
respectively, at the position of minimal mixing. Due to its MSSM nature the lightest Higgs boson 
receives large one-loop corrections. The correction is ~ 35 GeV for Mhi at its maximum value, 
which is taken at Mjj± ~ 515 GeV, and can even triple the tree-level mass at the borders of the 
Mff± range. The latter is mostly the effect of an already very small tree-level mass of ~ 30 GeV. 
The singlet-like H2 on the other hand receives smaller corrections which can nevertheless reach 10 
GeV. The investigated parameter range has been partially excluded as indicated in Fig. [9] (right). 
For Mjj± < 452 GeV light Higgs boson masses Mhi ^ 112 GeV have been excluded by the LEP 

^^See also Ref. [41] for a discussion of the Higgs boson masses and their dependence on the NMSSM parameters. 

'^^It should be noted that this is only approximately true for the matrix obtained via the procedure described in 
Sect. 12.51 In general this matrix, in contrast to the one obtained using — 0, is not unitary. As already mentioned 
the difference between the two approaches is small for our scenarios. 
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Figure 9: Left: The matrix element squared (7?.^)^ {i = 1,2) as function of A//f± at tree-level (dashed) and 
one- loop (full) for Hi (blue/dark grey) and H2 (red/light grey). Right: The corresponding tree- level and 
one-loop corrected masses with the exclusion limits. 



searches in ZH — )• Zbb [53]. The right exclusion limit Mjj± > 585 GeV is also due to the LEP 
exclusion of Mhi ^ 113.5 GeV in the ZH — )• Zbb channel. The Tevatron results do not exclude 
Higgs mass values. The exclusion limits obtained with HiggsBounds agree with ours. The maximum 
mass value taken by the MSSM-like lightest scalar Higgs boson is 140.5 GeV. It is just below the 
value excluded by the newest LHC limits [56]. In view of the uncertainties associated with the 1- 
loop Higgs mass corrections and also taking into account the fact that our exclusion limits are only 
a rough estimate, the scenario might be excluded for charged Higgs mass values around 515 GeV. 

The influence of the renormalisation scheme on the one- loop corrected Hi^2 masses is shown in 




Figure 10: The mass Mh^ of the lightest (left) and of the ncxt-to-lightest (right) CP-even Higgs boson 
as function of at tree- level (yellow/full) and at one- loop level adopting a DR (blue/dotted), a mixed 

(red/dashed) and an on-shell (green/small dotted) renormalisation scheme. 
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Fig. [ini They are shown as functions of Mh± in the OS, mixed and pure DR scheme and compared 
to the tree-level result. For the chosen parameter set the differences in the mixed and OS scheme 
are negligible in contrast to the DR scheme. This is to be expected as both in the mixed and 
the OS scheme Ax, respectively Mjj±, are renormalised on-shell, whereas in the DR scheme only 
the divergent part is included in the counterterm. The differences in the renormalisation schemes 
are, however, of maximally 0(10%) for the masses of the light Higgs boson. For the heavier scalar 
particles they are only ~ 1% of the one-loop corrected mass. 

In Figure [11] we show the tree-level and one-loop corrected masses of the pseudoscalar Higgs 
bosons as functions of Mjj± . Whereas for the heavier CP-odd boson the corrections are negligible 
they can be of 0(5%) for the lighter one at large values of Mjj±. Nevertheless they are much 
smaller than the corrections for the lightest scalar Higgs boson. 




425 450 475 500 525 550 575 600 425 450 475 500 525 550 575 600 

Mjj± [GeV] Mjj± [GeV] 



Figure 11: The mass Ma^ of the lighter (left) and of the heavier (right) pseudoscalar Higgs boson as a 
function of Mh± at tree- level (dashed/blue) and at one- loop (full/red). 



4 Summary and Conclusions 

In summary, we have calculated the one-loop corrections to the NMSSM Higgs masses in a renor- 
malisation scheme which mixes on-shell and DR conditions, by applying the latter solely to the 
parameters tan/3, v^. A, ^4^. We have compared our result to a pure DR and a pure OS-scheme. 
Apart from special parameter regions which unphysically blow up the counterterms, the results 
differ by at most 10% for the masses of the lightest scalar Higgs boson. Another estimate of the 
effect of the missing higher-order corrections is given by the variation of the renormalisation scale. 
We found the effect on the Higgs mass corrections to be 0(10%) or less. Altogether, the residual 
theoretical uncertainty at one-loop level can be estimated to be of the order of 10%. 

The bulk of the one-loop corrections stems from the top quark sector as it is known from the 
MSSM. This is reflected in the difference of the relative corrections for the lightest Higgs boson 
mass. The relative correction in the scenario, which we investigated, is equal to 55% when adopting 
the running DR mass and 79% in case of the top quark pole mass. Furthermore, the higher-order 



27 



corrections can shift the point of cross-over between singlet-Uke and MSSM-hke behaviour of the 
Higgs bosons. This dictates also the amount of couphng to the gauge bosons and fermions and hence 
influences the Higgs boson phenomenology. The precise knowledge of the higher-order corrections 
is therefore indispensable for a proper investigation of specific production and decay scenarios and 
a proper interpretation of the experimental results. 



Appendix 

A The CP-even and CP-odd Mass Matrices 

In the following we display the mass matrices in terms of the parameters for which we apply our 
renormalisation conditions. We repeat them here for completeness, cf. also Eq. ()29p . 

th^, th^, th,, e, M^, M|, tan/3, M|±, A, k, Vs, . (104) 

Note, that in the mass matrices we paid attention to keep the distinction in the angle (3 and the 
angle /3b- The difference of these two angles is denoted as 

A/3 = /3 - /3b . (105) 

Here, /3 is defined by the ratio of the vacuum expectation values of the neutral components of the 
two Higgs doublets Hu and H^, tan/3 = Vu/vd- The angle /3b on the other hand performs the 
rotation from the basis (0^,0^1, a^) to the basis {a,as,G) to separate a massless Goldstone boson. 
It is the angle /3 which receives a counterterm. The angle 13b also enters the scalar mass matrix. 
This is due to the replacement Eq. (|27p of Ax in terms of the parameter set Eq. (jl04p . 



The scalar 3x3 mass matrix in the basis = {h^, hu, hg)^ is given by the entries M. 



M|,, (i, j = 1,2,3), with 



2 2 

+ ^^-^ (106) 

,^2 gC/^Cfe , ,2 , . 1 £^rA#2 ,rA/r2 A/f2 v2 1 , A^M^S^g2/3 , . 

^Si2 - . ^2 i^hjptfjg + thj - --^[Mh± + [Mz - My^)c^p\ + -2 (107) 

^ A/3 '^^s'^Ap 

^ W W p P ^-^Qg^ 



+ ^^'^t'^'' (109) 
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■5 A/3 ^'^s'-^ij ^ 



e^Vs 



(110) 



Us n UgLj^p 

I -^W"^W"^2/3 r„ \ 2 7, ^-2 2 \ 2 2 n^-2 2 1 /i i i \ 

H -^^^ [2A M^rSy^S2l3 " kAc V,, - M^yB S2/3J . (Ill) 

e 

The entries Mp^. = ^p^^ (^) J = 1) 2, 3) of the pseudoscalar 3x3 mass matrix Mp in the basis 
= (a, as, G)'^ read 

n\2A^2 2 2 

M^,, = ""^^^"^ + M^, - M^ci, (112) 

2 

2 _ Myfiyf2^r,.2 ,^2,2 1 ^^^[+ + +2 1 



+ ^^^^^^^^[2AM2,.2^.2, - S^e'v's] (113) 

6 

Vz Us eUgC^i^ eUgC^i^ 

+ ^^^^^^[2\M^s'^s,p + ?>Ke\l] (115) 
e fg 

zeVsCj\i3 VsC/\p 

+ ^^^^P^[2XM^y,sl,s,, - ?,.e\l] (116) 



e^Vs 

f2 ^;„2 



j\,r2 ji#2 J. 2 AO I -^M^ ro\2 2 

Mp33 = M^± tan A/^ H — 2 [2A 



„2n 



e 



+ . ^2 [^fcdC/3-2/3B - ^h„g/3-2/3fl] • (117) 



B Fermionic Self-Energies 

In this Appendix we give the renormahsed fermionic self-energies which are needed in the determi- 
nation of the counterterms 5vs and 6k from the chargino and the neutrahno sector, respectively. 
We start from the general structure of a fermionic self-energy given by Eq. (|60p , which also applies 
to the corresponding renormahsed self-energies. The left- and right-chiral chargino fields V^lr ™ 
the interaction basis are given in terms of the Weyl spinors for the gaugino fields Wi , W2 and the 
charged components of the higgsino fields , , 
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with W"^ = [Wi =F iW2)/V2. The char gino fields in the interaction basis are replaced by the 
renormalised fields and the corresponding field renormalisation constants, 



n,R^[^+2^ZL,R)i^lR^ ^itl^ ^ZL,R=y Q^'"^ 5Zl,r,]' ^^^^^ 
and the 2x2 chargino mass matrix by 

Mc ^ Mc + 5{Mc) . (120) 

Performing the rotation of V'l to the fields Xl r the mass eigenbasis with the unitary matrices 
U and V, 

xt = V^t > Xn = U^R (121) 

yields for the renormalised self-energies T, in the mass eigenbasis in terms of the unrenormalised 
self-energies, field renormalisation constants and mass matrix counterterm 

Sj±(fe2) = J:^i{k^) + ^U*i5Zn + 5Z%)U^ (122) 

t^±{k^) = ^L^(k^) + ^V{6ZL + 6Zl)V^ (123) 

t^Uk^) = J:^i{k^)-^U*{5ZRMc + Mc6ZL)V^ -U*5{Mc)V^ (124) 

t^lik'^) = $]%(A:2) - IvidZlMl + m1.5Z*r)U^ - V5{mI.)U^ . (125) 

A 2 

The self-energies are 2x2 matrices. The corresponding renormalised self-energies in the neutralino 
sector, which are 5x5 matrices, are obtained analogously and read 

sJo(fe2) = ^^o{k^) + ^Ar*{szR + sz*j,)Ar^ (126) 

t^oik^) = J:^,{k^) + ^Mi6ZL + 5Zl)M^ (127) 

t^S{k^) = J:^^{k^)-^M*{6ZRMN + MN6ZL)M^ -J\f*6{MNW^ (128) 

Sjo^(A:2) = ^^S{k^)-^M{6ZlMl, + Ml6Z*R)J\f'' -M6{Ml)J\f'' , (129) 

where is the unitary 5x5 matrix which performs the rotation from the interaction basis to 
the neutralino mass eigenbasis. Allowing also for negative neutralino mass values, as we do, the 
matrix is real. The neutralino mass matrix Mjv is given in Eq. (jl9p and SZl^r denote diagonal 
5x5 matrices with the field renormalisation constants in the interaction eigenbasi^^, 5Zi^^r. 
{i = 1, 5), as entries. Note that 6Zl and 6Zr are related due to the Majorana character of the 
neutralinos. 



^^For better readability we did not choose a different notation for the field renormalisation constants in the chargino 
and the neutralino case. They are understood to be different, though. 
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